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Abstract. In the present survey we present some of the recent results 
concerning the geometry of nodal lines of random Gaussian eigenfunc- 
tions (in case of spectral degeneracies) or wavepackets and related issues. 
The most fundamental example, where the spectral degeneracy allows 
us to consider random eigenfunctions (i.e. endow the eigenspace with 
Gaussian probability measure), is the sphere, and the corresponding 
eigenspaces are the spaces of spherical harmonics; this model is the pri- 
mary focus of the present survey. The list of results presented is, by no 
means, complete. 

1. Introduction 

Nodal patterns (first described by Ernest Chladni in 18th century) ap- 
pear in many problems in engineering, physics and the natural sciences: they 
describe the sets that remain stationary during vibrations, hence their im- 
portance in such diverse areas as musical instruments industry, mechanical 
structures, earthquake study and other areas. They also arise in the study 
of wave propagation, and in astrophysics; this is a very active and rapidly 
developing research area. 

1.1. Some basic notation. Let {M,g) be a compact Riemannian surface 
(for example S^, the two dimensional unit sphere equipped with the round 
metric), and A be the Laplace-Beltrami operator on M. We are interested 
in the eigenvalues A and the corresponding eigenfunctions (p of —A, so that 

A4> + X(j) = 0. 

In case Ai has a boundary, we impose either the Dirichlet boundary condi- 
tion 4>\dM = 0, or the Neumann boundary condition 

^"^1 -n 
ou 

or any mixture of the conditions above. The general spectral theory states 
that there is a complete orthonormal basis of {M ) which consists of eigen- 
functions, i.e. we may choose a sequence of functions 

and corresponding nondecreasing sequence of eigenvalues {Xj}°^i so that 
the orthonormal set {(pj} spans the whole of L^(A^). Note that we allow 
multiple eigenvalues i.e. spectral degeneracies. 

Let (/> : 7W — 7- M be any real valued function. The nodal line of is its 
zero set 

(/.-i(O) = {x £M: (pix) = 0}. 
1 
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In principle, a nodal line might have self-intersection^; however generically 
it is a smooth 1-dimensional curve with components homeomorphic to either 
the circle ("closed component") or an interval (in this case it must intersect 
the boundary; it is called an "open component"). We are interested in the 
geometry of the nodal lines of (pj as j — )• oo. The most basic aspect of the 
nodal line is, of course, its length; let us denote Ij to be the length of the 
nodal line of (pj. In this survey we will only consider the 2-dimensional case; 
however most of the results presented extend for higher dimensions. 

1.2. Yau's conjecture and Berry's RWM. Yau conjectured [231 [23] that 
for any smooth M, Ij are commensurable to y^Aj for any smooth metric g 
in the sense that there exist two constants c{Ai), C{M.) so that 

(1) ciM)^/^] <lj <C{M)^/X'j 

for every j > 1- The lower bound was proved by Bruning and Gromes [9] and 
Bruning [8] for the planar case. Donnelly and Fefferman [12] finally settled 
Yau's conjecture for real analytic metrics. However, in its full generality, 
Yau's conjecture is still open. 

In his seminal work [3j , Berry argued that the high energy behaviour of the 
eigenfunctions should be universal, at least for "generic" chaotic surfaces M. 
He proposed to compare an eigenfunction with eigenvalue A to a "typical" 
instance of an isotropic, monochromatic random wave with wavenumber 

A; = 

(nowadays called Berry's Random Wave Model - RWM). A 1-dimensional 
version of the random wave was used by Rice in order to investigate the 
likelihood of a given signal to exceed a level. Longuet-Higgins generalized 
Rice's model to 2-dimensional plane to describe the movement of the sea 
and ocean waves. 

There are several ways to construct the ensemble of random waves as- 
sociated with wavenumber k. One way to do it is consider summations of 
typeO 

on M?, where 9j are random directions drawn uniformly on the unit circle, 
and (pj G [0, 27r) are the random phases. One would like to define the random 
wave Uk{x) on as the limiting ensemble 

Uk{x) = lim Uk-j{x); 

which should converge in distribution. Another mathematically rigorous 
way to define the wavenumber k Random Wave is to identify it as the 
unique Gaussian isotropic random field (ensemble of functions) with co- 
variance function 

(2) rRWM{x,y) = Mk\x -y\), 



For example the eigenfunction (f>{x,y) — sin(107rx) sin(207ry), defined on tlie torus 
= m} having eigenvalue A = SOOvr^. 
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where Jq is the usual Bessel function. 

Since, according to the RWM, the random waves model the high-energy 
eigenvalues, the nodal lines of random waves should also model the nodal 
lines of honest eigenfunctions. This approach allows us to study local quan- 
tities like the nodal length, boundary intersections and intersections with a 
test curve etc. 

Suppose, for example, we are interested in the nodal length on the torus. 
Then we are to choose a representative planar domain U <^ M? (e.g. a 
rectangle with the same aspect ratio and area as the torus), and study the 
distribution of Zjj.^, the nodal length of random wave with wavenumber 

k = \/A inside U. It is easy to compute the expected length to be of order 
of magnitude 

E[Zj^.^] ~ const -Vxpl, 

where \U\ is the area of \U\, and Berry argued [¥] that the variance should 
be of order 

(3) Var (2:^.^) ~ const ■ \U\ log A. 

1.3. Bogomolny and Schmit's percolation model. The RWM, how- 
ever, does not help if one is interested in making predictions regarding the 
more subtle (and arguably, more interesting) aspects of the nodal structures 
such as the number of nodal domain^, their size distribution, the size of 
the largest nodal domain, the inner radius etc. For this purpose an elegant 
independent bond percolation-like model was introduced by Bogomolny and 
Schmit According to this model, the nodal domains should correspond 
to the clusters connected by open bonds, and the nodal line corresponds to 
these clusters' boundaries [7]. Let Uj = N{(j)j) be the number of nodal do- 
mains of and N{u^; U) the number of nodal domains of random wave 
u^.jj on U. The main criticism against this model is that the independence 
assumption ignores all the dependencies that occur between the bonds; to 
try to justify the independence assumption, Bogomolny and Schmit [6] apply 
a heuristic principle, the so-called Harris criterion. 

By the classical Courant Nodal Domain Theorem (see e.g. |10]). vj < j, 
and Pleijel [l19J asymptotically improved the latter to 

limsup^ < 0.691.... 

j— s>oo J 

On the other hand, no nontrivial lower bound for Vj could be found, since 
one may find a sequence of energy levels Aj^. on the torus (say), so that the 
corresponding eigenfunctions would have only 2 nodal domains. Bogomolny 
and Schmit [5] used the general percolation theory to predict that 

iV(n^; U) 



In reality the summation is slightly more complicated than the one presented - see 
e.g. 0, (1). 

^The nodal domains are the connected components of the complement of the nodal 
line. 
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should be asymptotically Gaussian, with mean and variance proportional to 

\u\ ■ [vxy = \u\ ■ X. 

More strikingly, in their later paper, Bogomolny et al. [7j argue that since, 
according to the recent developments in the percolation theory (see e.g. 
Smirnov [21]), the "interface" (cluster boundaries) should converge to SLEq, 
the distribution of the largesl|jcomponent of the nodal line (rather than its 
length) should converge to SLEq. 

1.4. Equidistribution conjecture. It is conjecture41 that on any chaotic 
surface M (for example, any negatively curved surface, or any ergodic bil- 
liard), the nodal lines are asymptotically equidistributed in A^, so that, in 
particular, the nodal length has asymptotic shape of type 

(4) ij ~ CM ■ yA~ 

(refinement of Yau ([T])), for some constant cm > 0. Despite the fact that 
heuristically, ^ follows from the RWM (it follows for example from the 
results mentioned in the end of Section [TT2|), this conjecture seems extremely 
difficult or even out of reach by the present analytic methods, and it seems 
highly unlikely that it is going to be settled in the near future. 

However, some information could be inferred from the completely inte- 
grable case, even though the picture that emerges here is very different. For 
instance, one may use the spectrum degeneracy of the standard torus to eas- 
ily construct sequences of eigenfunctions (pm-j and (/)„2.j 5^ = 1,2,^ = 1,2,..., 
so that 

length (</'4(0)) ^Ci-^^, 

i = 1,2, with ci 7^ C2; one may obtain such sequences on the sphere only 
slightly modifying the same argument. One way to infer some information 
is use the following heuristic principle: 

Principle 1.1 ("Word exchangeability"). Any property satisfied by generic 
eigenfunctions on (all) completely integrable manifolds is also satisfied by all 
eigenfunctions on a generic chaotic manifold. 

1.5. Acknowledgements. The author would like to thank Zeev Rudnick 
for suggesting to consider some of the problems in this survey, and having 
many deep and fruitful conversations. I would also like to thank Mikhail 
Sodin for many extremely stimulating and fruitful discussions regarding 
some of the subjects that appear in this survey and their context, and many 
useful comments on an earlier version of this manuscript. In addition, it 
is important for me to acknowledge the organizers of the Dartmouth Inter- 
national Conference in Spectral Geometry for organizing such a wonderful 
conference, that included many wonderful talks and speakers, in addition 
to the free and relaxed informal environment or atmosphere, stimulating 
and encouraging collaboration and new ideas' exchange, and also for the 
generous financial support. 

^According to the percolation theory, there exists exactly one component that "covers" 
the whole domain; this is the only macroscopic component of the nodal line. 

^In a recent survey by S. Nonnenmacher ,_17„, this conjecture was attributed to S. 
Zelditch. 



on the nodal lines of random and deterministic laplace eigenfunctions 

2. Some results 

2.1. Spherical harmonics. It is well known that the eigenvalues E of the 
Laplace equation 

Af + Ef = 

on the 2-dimensional sphere 5^ are all the numbers of the form 

(5) En = n{n + 1), 

where n is an integer. The corresponding eigenspace is the space Sn of 
spherical harmonics of degree n; its dimension is 

J\fn = 2n+ 1. 

Given an integer n, we fix an L?'{S'^) orthonormal basis of £n 

r/^(x),...,r/^„+i(x), 

giving an identification £n — K-^". For further reading on the spherical 
harmonics we refer the reader to [1], chapter 9. 

2.2. Random models. In case of spectral degeneracy, such as the sphere or 
the torus, we may consider a random eigenfunction lying inside an eigenspace. 
For the sphere, we define it as 

(6) fn{x) = \Jj^Yl "'kVkix), 

where ak are standard Gaussian A^(0, 1) i.i.d. That is, we use the identifi- 
cation 

to endow the space £n with Gaussian probability measure v as 

where a = (oj) G E/^" are as in ([6]). 

Note that v is invariant with respect to the orthonormal basis for Sn- 
Moreover, the Gaussian random field /™ is isotropic in the sense that for 
every xi, . . .xi G 5^ and every orthogonal R S 0(3), 

(7) ifniRxi), . . . , fniRxi)) = (/„(xi), . . . , fnixi)) . 

There exists yet another way to define /„: it is the centered Gaussian 
isotropic random field with covariance function 

(8) rn{x,y) := E[/„(x) • fn{y)] = Pn{cosd{x,y)), 

where P„ are the well-known Legendre polynomial of degree n, and d is the 
spherical distance. The Legendre polynomials admit Hilb's asymptotics 



(9) p^(cos((/<)) « W-4Jo(0(n + l/2)), 



sm ( 



i.e. almost identical to RWM ([2]), up to the "correction factor" y ^j^- This 
factor seems to "keep a trace" or "remember" about the geometry of the 
sphere. 
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For generic manifolds there are no spectral degeneracies, so that we 
should consider linear Gaussian combinations of individual eigenfunctions 
( "wavepackets" ) . The two most accepted models are the so-called long en- 
ergy window and short energy window. In the long window case we consider 
Gaussian combinations of eigenfunctions with eigenvalue lying in the window 
[0, A] for A — ;> oo 



where x G A^, and 

Nm{\) = *{\j<\} 

is the spectral function; the reason we took the square root of the eigenvalues 
is that it makes it more convenient to write the short energy window random 
function 

/f (x) = , ^ = aj(j)j{x). 

jNMii^+m-NMW •/X<^/x]<•/X+l 

The short energy modeQ is considered more significant, as it is more repre- 
sentative of the individual eigenfunctions Q; however working with the long 
energy window is relatively easier. The spectral prefactor in the definition 
of f^'^ was introduced to make the expected L^-norm unity. 

As usual, for any random variable X, we denote its expectation EX. For 
example, with the normalization factor in 1^, for every n fixed point x G 5^, 
one has 

(10) E[/„(x)2] = ^j;r?,"(x)2 = l, 

k=i 

a simple corollary from the Addition Theorem (see [I]). 
Any characteristic X(L) of the nodal line 

f-HO) = {x G : /„(x)=0} 

is a random variable. The most natural characteristic of the nodal line of 
/„ is, of course, its length Z{fn)- One may then study the distribution of 
the random variable Z{fn) for a random Gaussian G En, as n ^ oo. It is 
also natural to consider the number N{fn) of the nodal domains of /„, i.e. 
the connected components of 

its distribution should be consistent to the one predicted by Bogomolny and 
Schmit based on their percolation model (see Section [1.31 above) . 



^The window [A, A + 1] may be replaced by [A, X + a] for any constant a > 0. 

n+Vn 

Much like the random trigonometric polynomials Yl o-n cos{nt) on [0, 2n] is more 

n=JV 

JV 

representative of cos{Nt) than a„cos{nt), where in both summations a„ are standard 

n=l 

Gaussian i.i.d.; for example it possesses asymptotically the same number of zeros. 
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2.3. Some generic results. Berard t2j, and subsequently Zelditch [28] 
found that the expected nodal length of the long energy window random 
functions is 

E [^(/f)] ~ const ■ VX, 

consistent with Yau. Zelditch also extended this result to the short 
energy window case 

E )] ~ const ■ \/A. 

In addition, for g real analytic, Zelditch |28j considered the complexified 

manifold (A^C)5c) (whose projection on $52 = is {M,g)), the analytic 

(P C"-L S 

continuations (pj, and the corresponding random combinations fx' ' , de- 

fined analogously to the real random combinations ' . In this case the 
zeros are isolated points in A^ci Zelditch [28j found that their expected 
number is again proportional to \/A- 

Toth and Wigman [22] considered the number of boundary intersections 



L S 

of the nodal line of ' in case is a generic billiard (i.e. a planar oval 
with a smooth boundary), or equivalently (up to the factor 2), the number 
of open components. They found the correct order of magnitude for the 
expected number of intersections to be 



^(/a'^)1 ~ const\dM\V\, 



E 

where \dM\ is the boundary length of 7W, and the constants differ in the long 
and short window cases. This result is consistent to both Yau's conjecture, 
random wave model, and the resulting interpretation of the boundary trace 
as approximating trigonometric polynomials; the asymptotics depends only 
on the boundary length of the billiard, notably independent of its shape. 

2.4. Number of nodal domains. Nazarov-Sodin [18] found the correct 
order of magnitude for the expected number of nodal domains of random 
spherical harmonics, and established an exponential decay result for devia- 
tions from the mean. 

Theorem 2.1 (Nazarov-Sodin |18]). There exists a constant a > so that 
the expected number of nodal domains is asymptotic to 

(11) E[N{fn)] = an^ + o{n^). 

Moreover, for every e > 0, there exist two constants c(e),C(e) > 0, so that 

Nifn) 



(12) Pr 



> e < C{e)e 



-c(e)n 



The result pip on the expected number of nodal domains is of more 
general nature: it extends to a wide range of sequences of random fields 
(to appear in a paper by Nazarov and Sodin). For example, rather than 
taking a random element lying in a single spherical harmonics space, one 
may superpose elements from several spaces). However, unlike the rapid 
decay (jl2p in the particular case of spherical harmonics, in the more general 
situation, Nazarov-Sodin's result does not prescribe the rate of decay of the 
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tails of the distribution. Instead, they prove the weaker statement: for any 
e > 

Nifn) 



Km Pr 

n— >oo 



> e = 0. 



One disadvantage of the results above and the method of their proofs is 
the fact that the constant a > 0, whose existence is established, remains 
mysterious and completely open; one cannot establish the dependence of a 
depends on the underlying random field. As an example, one may oppose 
fn to Qn, a superposition of random spherical harmonics of degree < n. A 
generalized version of Nazarov-Sodin's Theorem implies the existence of a 
constant so that 

Theorem 12.11 and its generalization do not shed a light on the the relation 
between a and b. 

Though the expectation result is consistent with the percolation model, 
Theorem 12.11 gives us no clue what would be the result for the variance. 
However, it provides us with a very strong rate of decay (namely exponen- 
tial); the decay results are usually complementary to the variance results. 
The authors also proved that the prescribed rate of decay cannot be im- 
proved, so that the exponential decay they establish, is of the correct order 
of magnitude. 

2.5. Nodal length of random spherical harmonics. It is a standard 
application of the Kac-Rice formula (see e.g. [IT]) to compute the expected 
nodal length of random spherical harmonics [2] 

(13) E [Z{fn)] = c • v^, 
where 

c = V^vr, 

(see also [TB] and [25]). Our main concern in this pursue is the subtle 
question of the variance. 

Based on the natural scaling of the sphere (e.g. the relation to the Le- 
gendre polynomials, in particular ([8|)), we conjectured [25j, that 

Var(Z(/n)) ~ const • n. 

Surprisingly, the variance turned out to be much smaller, due to an unex- 
pected cancellation ("Berry's cancellation phenomenon"). We derived the 
following asymptotics for the nodal length variance, improving the earlier 
bounds of Neuheisel [16] and Wigman [25j: 

Theorem 2.2 (Wigman [26j). As n ^ co, one has 

(14) Var(^(/„)) = -logn + 0(l). 

Note that the leading constant || in (|14p is differenlQfrom the one pre- 
dicted by Berry for the RWM (see ([3]) ) . Our explanation for this discrepancy 
is the nontrivial local geometry of the sphere. It seems reasonable that for 

L S 

a generic chaotic surface, the nodal length variance for ' should be log- 
arithmic; the leading constant is then an artifact of the local geometry. 
One of our central goals is to find this dependency explicitly, namely, given 
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a Riemannian surface, {Ai,g), compute the leading constant in front the 
logarithm (if this prediction is indeed correct). 
Theorem 12.21 implies that the series 



of variances of the normalized length 

2ifn) 



nzifn 



is convergent. Together with the Borel-Cantelli Lemma, it implies that for 
independently chosen /„, 



En 



almost surely, where c > is the same constant as in ()13p . 

The same problem may be also considered on the standard 2-dimensional 
torus T = M^/Z^. Here the eigenvalues are of the form 

£'J = 47r^n, 

where n is an integer, expressible as a sum of two integer squares, and 
the corresponding eigenspace is spanneclBby functions cos {2tt{X, x)) and 
sin (27r(A, x)), where A G with ||A|p = n, are all the lattice points lying 
on the circle of radius y/n; its dimension is r2{n), the number of represen- 
tations of n as a sum of two squares. The Gaussian random eigenfunction 
is a stationary random field with the covariance function 

(15) rl{x) = ^ Yl cos(27r(A, x)). 

l|A|P=n 

It is again standard to compute that the expected nodal length of this 
ensemble to be proportional to \/eX \2Qh we are interested in the 
asymptotic behaviour of the variance again. This question was initially con- 
sidered by Rudnick and Wigman [20]; however it got only a partial answer 
then. An (almost) complete answer will be given in the forthcoming paper 
Krishnapur-Kurlberg- Wigman [13j . 

Even though we have an explicit expression (jlSp for the covariance func- 
tion, no analogue of ^ is known for the asymptotic long-range behaviour 
of (recall dH])). As a replacement, we cope with some subtle issues of 
the arithmetics of lattice points lying on a circle. Here as well we observed 
the "arithmetic Berry's cancellation" , a phenomenon of different appearance 
but similar nature to "Berry's cancellation phenomenon". 



Since /„ is odd for odd n and even for even n, the nodal lines are invariant w.r.t. the 
involution x i— >■ —x. Therefore the natural planar domain to compare would be one of area 
of a hemisphere rather than of the full sphere. 

^Note the invariance w.r.t. A i— > —A, so that we need to factor the set of lattice points 
by ±. 
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2.6. Level exceeding. Let us define tlie spherical harmonics level exceed- 
ing measure as follows: for all z £ (—00,00), 



(16) ^niz):= / lif nix) < Z)dx, 

where l(-) is, as usual, the indicator function which takes value one if the 
condition in the argument is satisfied, zero otherwise. In words, the function 
^niz) provides the (random) measure of the set where the eigenfunction lie 
below the value z. For example, the value of ^n{z) at 2; = is related to the 
so-called defect 

Vn := meas (/~^(0,oo)) - meas (/~^(-oo,0)) 

by the straightforward transformation 

Vn = 4TT- 2«>„(0). 

Of course, 47r — ^n{z) provides the area of the excursion set 

An{z) := {x : fn{x) > z} . 

Clearly, for all z eR, 

E[$„(z)] =47r$(z), 

where $(•) is the cumulative distribution function of the standard Gaussian. 
The following lemma (see Marinucci-Wigman [Ij]) deals with the variance 
of ^n{z) as n — )• 00. 



Lemma 2.3. For every z G 



Var($„(z)) = z^</.(z)^--+0, , , 

n \ 



where (p is the standard Gaussian probability density function. 

In particular, for z ^ 0, Lemma 12.31 gives the asymptotic form of the 
variance as n — )■ 00. In contrast, for z = (this case corresponds to the 
defect), this yields only a "o" -bound and one needs to work harder to obtain 
a precise estimate; we do so in the forthcoming paper: 

Theorem 2.4 (Marinucci-Wigman [15j). As n ^ 00, the defect variance is 
asymptotic to 



Var(P„,)~C-^, 



where C > is some constant. 



In light of Lemma 12.31 it is then natural to normalize $n(-2) and define 
the spherical harmonics empirical process by 



(17) Gniz) := 



[ 1 Unix) <z)dx- {47r X $(z)} 

^52 



for n = 1, 2, z G (—00, 00). 

In |14] we proved the following result: 
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Theorem 2.5 (Marinucci-Wigman [14]). (The Uniform Central Limit The- 
orem) As n ^ oo, the process G„(z) converges in distribution to Goo{z), 
where Goo (z) is the mean zero, degenerate Gaussian process on M given by 

Goo{z) = Z(j){z)Z 

with Z ~ iV(0, 1) standard Gaussian random variable. 

This resuh, in particular, imphes the fuh asymptotic dependence of Gn{z) 
for different values of z, as n — )• oo. See the next Section for some explanation 
to this phenomenon. 

2.7. Nodal line vs. Level curves. Interestingly, the behaviour of level 
curves f^^{L) for L 7^ is very different compared to the behaviour of nodal 
lines. Let Z^{fn) be the length of the level curve f~^{L). It is standard to 
compute the expected length, using the Kac-Rice formula 

2 , 



consistent with the nodal case L = 0. However, unlike the nodal lines, level 
length variance is asymptotic to [27] 



(18) Var(^^ (/„)) ~ caL^e-^ • n; 

it is also interesting to observe the fact that the leading term depends on 
-L^ (a priori, the dependence on L should be symmetric w.r.t. L 1— t- — L, 
however we would rather expect L^; its dependence cancels out - another 
obscure cancellation related to this problem). 

Moreover, the length of the level curves is asymptotically fully correlated, 
in the sense thalE] 

(19) p(Z^l(/n),^^'(/n)) = 1 - On^oo (1) • 

Let us relate between the latter and the setup of Theorem 12.51 One may 
express ^n{z) in terms of the level lengths using 

z 

(20) ^n{z)= j Z''{U)dL. 

— oo 

Intuitively, the asymptotic degeneracy of ^n{z) for different values of z is 
then an artifact of the asymptotic full dependence ()19p of the individual 
values in the integrand on the RHS of (pOj) . 

One possible explanation for the phenomenon ()19p is the following conjec- 
ture, due to Mikhail Sodin (see Marinucci-Wigman [15] for further reading 
on this conjecture). For x G 5^ and L G M let 

■^x — ■^x ifn) 

(the "local length") be the (random) length of the unique comp onenlF%f 
f~^{L) that contains x inside (or 0, if /„ does not cross the level L). 
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For two random variables X, Y the correlation is defined as p{X, Y) 
Cov(X Y) I I ^ ^ measures the linear correlation between X and Y. 

^^We should assume that fn{x) ^ L; the latter is satisfied almost surely. 
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Conjecture 2.6 (M. Sodin). The local lengths are asymptotically fully de- 
pendent in the sense that for every x £ S'^ and -^1,^2 £ 

Intuitively, it should be clear that Conjecture 12.61 implies (fTUj) (and thus 
also the asymptotic degeneracy of the level exceeding measure via (pO]) ). 
since is some summation of Z^ over some x on the sphere. 

All the results above piece up nicely together while dealing with the cor- 
responding questions in case of Gaussian subordinated random field£3: let 
G : M — 7- R be a (possibly nonlinear) nice function and define the random 
field gn by 

gn{x) = G{fn{x))] 

gn is Gaussian subordinated. We are interested in the nodal length of gn 

Zg„ =length(5-i(0)), 
so that if zi, . . . , Zk are all the zero S)f G. It is obvious that 

k 
i=l 

Therefore, the expected nodal length of gn is 

k 
i=l 

for some explicitly given ci > 0. 

It is not difficult to see that (fTSj) together with (fT9]) gives an elegant and 
compact asymptotic result (for n — )• 00) for the nodal length variance of gn 
as 

/ k \ 2 



vi=l 



with some explicit C2 > 0, provided that Zi ^ for at least one index i. 
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